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Abstract We present an overview of recent developments, based on on-shell techniques, in the 
calculation of multi-parton scattering amplitudes at one loop that are relevant for phenomeno- 
. logical studies at hadron colliders. These new on-shell methods make efficient use of the physical 

C^") ' properties of the hard scattering, such as unitarity and factorization. 
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1 Introduction 



With the recent first collisions at the LHC we are entering a new era of discovery 
in particle physics. Colliding protons at very high energies makes the LHC a 
fertile environment for the production of high-multiplicity events. If we are to 
take full advantage of the discovery potential at high energy hadron colliders such 
as the Tevatron and the LHC, we need to have a precise understanding of the 
physics that will occur there. This necessitates computations of the Standard 
Model background, especially of QCD processes, to at least next-to-leading order 
(NLO) in the perturbative series. 

Historically, the bottleneck in NLO computations has been the one-loop vir- 
tual contributions. Over the last few years rapid progress has been made in the 
development of new techniques for these one-loop computations. These advance- 
ments have been motivated both by a greater desire to understand the structure 
of scattering amplitudes as well as the need for improved efficiency and automa- 
tion in the computation of one-loop matrix elements. These matrix elements are 
needed, for example, for the precise computation of many background processes 
at the LHC (see e.g. ([I])). The ability to automate and "mass-produce" ampli- 
tudes requires approaches that are both numerically stable and straightforward 
to implement as an algorithm. 

The standard approach to performing this class of computations has heavily 
relied upon Feynman diagram techniques. There have been many impressive 
results with this approach (see ([!]) and references therein; for example, cross sec- 
tions for 6-point processes that have been computed via Feynman diagrammatic 
methods include ([!])). However, Feynman diagrams suffer from two problems. 
For one, there is a factorial growth in the number of terms as the multiplicity of 
partons in the process increases. Furthermore, each Feynman diagram is gauge 
dependent. This means that there will be large cancellations between terms that 
combine to give the gauge-independent amplitudes. It is these two problems that 
make automated approaches using Feynman diagrams difficult as the number of 
partons increases. 

Consequentially, the focus of much recent progress has been to side-step these 
issues. On-shell recursion and unitarity methods work with gauge-independent 
amplitudes as building blocks instead of Feynman diagrams. It is these new 
techniques that we focus on in this review. An earlier review of these techniques 
was presented in ([3]). However, many of the details presented there have been 
superseded by even more efficient techniques which we present in the following. 

In general, as we will explain in more detail below, a one-loop amplitude can be 
decomposed into a set of scalar box, triangle, bubble, and tadpole integrals, that 
is, integrals with four, three, two, or one loop propagators, respectively. These 
scalar integrals contain all the logarithmic and polylogarithmic dependence of the 
amplitude and are multiplied by rational coefficients. In addition, there are also 
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purely rational terms in the amplitude. 

The original unitarity approach of applied two-particle cuts in four di- 
mensions and was used to produce many results ([5]). The terms containing 
(poly)logarithms and associated constants could all be computed via two-particle 
cuts. Terms with only a rational dependence on momentum invariants however 
could not be computed in this way and separate techniques were required ([!]). A 
systematic approach for the computation of these rational terms without the use 
of Feynman diagrams did not appear until a few years ago. The starting point 
for these developments were the recursion relations, developed at tree level by 
Britto et al. (BCFW) ([6j). An earlier version of a tree- level recursion relation 
was the off-shell Berends-Giele recursion ([7]). The proof of the BCFW on-shell 
recursion relations only relies upon the factorization properties of the amplitudes 
and on Cauchy's theorem. They could therefore be adapted to the more compli- 
cated problem of computing the rational parts of loop amplitudes (j8,9). On-shell 
recursion for the rational parts of loop amplitudes was first used in an analytic 
context (|9j) and then further adapted into a numerical procedure (fT0|) . 

At the same time, improvements to the original unitarity approach were also 
occurring. Brandhuber et al. used two-particle cuts with MHV vertices to com- 
pute certain sets of one-loop amplitudes (jlip . The application by Britto et al. 
of generalized unitarity ([5l ll2p to the computation of box coefficients highlighted 
the benefits of examining not the one-loop integral as a whole, but its integrand. 
The work of Ossola, Papadopoulos and Pittau (OPP) ([TBI [Til) followed in this 
vein. Upon separation of the integrand into a standard set of basis terms, the 
problem could be reduced to solving for their coefficients numerically. The simple 
analytic extraction of triangle and bubble coefficients was the focus of the work 
by one of the authors (15), where the known analytic behavior of the integrand 
was used to straightforwardly extract the bubble and triangle coefficients. A nu- 
merical adaptation of this procedure suitable for automation was then presented 
in Ref. ([TO]) . 

The investigation of the one-loop integral itself has also presented new direc- 
tions to explore. Britto et al. (|16p showed how the integral when written in a 
canonical form can be directly integrated. This procedure was developed fur- 
ther to produce new results ()17p and also to provide a general analytic structure 
for the coefficients of one- loop amplitudes (I18pi9p . Taking further advantage of 
the analytic properties of the two-particle cut amplitude, Mastrolia (|20p applied 
Stokes' theorem and a generalized residue theorem to compute bubble coefficients 
via direct integration. 

Extending the four-dimensional cut techniques to D dimensions (in dimen- 
sional regularization) has also provided a second very fruitful approach for the 
computation of the rational terms. In D dimensions the rational terms develop 
branch cuts and are so accessible via unitarity cuts. Giele et al. (|2ip elucidated 
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the extra structures present in D dimensions beyond those of the original four- 
dimensional OPP integrand. Hence they were able to compute the rational terms. 
Following up on this and an earlier result relating masses to D-dimensional uni- 
tarity (]22p . Badger presented an alternative computational approach where the 
(D — 4)-dimensional terms are treated as an additional mass in the loop (|23p . 
A different approach advanced by Papadopoulos et al. (]24[) involves splitting the 
rational computation into two pieces. One part is computed from the one-loop 
integrand and is an extension to the original OPP approach, while the second 
part comes from a reduced form of Feynman diagrams. 

These approaches have been implemented in several automated tools for the 
computation of one-loop amplitudes: BlackHat (|10p . CutTools/OneLOop (j25|126j) . 
Rocket (|27p . and others (28). These programs, combined with tools for the 
real emission part (|29U3Up have yielded a host of new results at next-to-leading 

order dHUESHHalEllESlES]) - 

Below we will review the main ideas of all these developments. Due to lack of 
space, we regret that we can neither present all details nor an exhaustive list of 
all results and refer the reader to the cited literature and references therein. 

We begin with a review of our notation and the general structure of multi- 
parton scattering amplitudes and explain how to construct such amplitudes re- 
cursively at tree level. We will use the spinor formalism, and those readers 
unfamiliar with spinor techniques and their application to the calculation of tree- 
level and one- loop amplitudes may wish to consult for example Refs. (137 p 38 1l39|) . 
Section [3] discusses the use of (generalized) unitarity to construct one- loop ampli- 
tudes from tree level amplitudes, both the cut and the rational part. In Section [5] 
we explain how to alternatively employ a recursive approach to obtain the purely 
rational non-logarithmic terms of amplitudes that cannot be constructed via uni- 
tarity in four dimensions. We conclude with a summary and give an outlook on 
the expected progress in the near future. 

2 Structure of Amplitudes 

Below we briefly review our notation, which closely follows Ref. (38) (see also 
(f3"7|) ). We then discuss the general structure of tree and one-loop amplitudes in 
renormalizable gauge theories to set the stage for the subsequent sections which 
discuss new methods for the computation of these amplitudes. 

2.1 Notation and Color Decomposition 

In the following we will consider amplitudes where all coloi0 and coupling infor- 
mation has been stripped off. We can express any amplitude in terms of some 
basic color (and coupling) factors which are multiplied by color-ordered subam- 

1 Color here and below refers to any group theory factors. 
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plitudes, or primitive amplitudes. These primitive amplitudes, which are defined 
for specific cyclic orderings of the external partons, carry all the kinematic in- 
formation but no explicit color indices. The full amplitude is then assembled 
from the primitive amplitudes by dressing them with appropriate color factors. 
There exist different, but equivalent, ways of grouping the primitive amplitudes 
together into so-called partial amplitudes, and we refer the reader to the litera- 
ture for further information, see (37, 38, 40l l4ip and references therein. In what 
follows below, we will discuss the computation of the primitive amplitudes, con- 
centrating on the description of the evaluation of the kinematic part. As shown 
recently in Ref. (|42p . the algorithms that we will present below can be extended 
to amplitudes with color information. 

We express the primitive amplitudes in terms of spinor inner products j§ 

(jl) = {r\l + )=^(k J )u + (k l ), [jl] = (j + \l-)=u + (k J )n.(k l ), 

(1) 

where u±(k) is a massless two-component (Weyl) spinor with momentum k and 
positive or negative chirality, respectively, which we also write as, 

(Aj) Q = [u + (ki)] a , (Xi)a = [u~(ki)] & . (2) 

Massless four-momenta can be reconstructed from the spinors by, 

k i( a li)cea = ($i) a a = (K)a(^i)a • (3) 

Spinor products can thus be used to construct the usual momentum dot products 
via 

(i j) \j *] = 2 Tr [titj] = 2k i ■ k j = s v ■ ( 4 ) 

Furthermore, we use the following notation for sums of cyclically-consecutive 
external momenta and their invariant masses, 

K., = # + + ■ + + (5) 

Si...j = At.,- (6) 



where all indices are to be understood mod n for n-particle amplitudes. 

The above formalism can be extended to include massive spinors and vectors, 
using the well-known decomposition of any, not necessarily light-like, four-vector 
k into a sum of two light-like four- vectors: 

Here q is a fixed light-like four-vector, fixing the axis of the spin for spinors, and 
k is the associated projection of the massive vector k. Massive spinors can be 



2 Note that we use the sign convention of most of the QCD literature, in the "twistor" 
literature a different sign convention for [j I] is used, for example in Refs. (|6|>. 
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constructed via 



and similarly for the conjugate spinors (|43p . Here, the label ± indicates that 
the spinors u± are eigenstates of the projector (1 ± ^7 5 ), with the spin vector 
s = k/m — m/{k ■ q) q. The normalization is chosen to allow a smooth limit to 
the massless case. 



2.2 On-Shell Recursions at Tree Level 

An efficient recursive technique for computing tree-level multi-parton scattering 
amplitudes was developed more than 20 years ago ([7]) and adapted for numerical 
implementation in various computer codes (|414I441[4"5|) . Berends-Giele recursion 
connects smaller off-shell currents together to produce amplitudes. More recently 
it was realized that through the use of complex kinematics amplitudes can be 
computed entirely using only smaller on-shell amplitudes. This leads to more 
compact analytic expressions not only at tree level but for rational terms also at 
loop level. Here we briefly review the on-shell recursion relations for tree level 
amplitudes found and proved in Refs. ([6|). A recursive approach at loop level is 
not quite so straightforward, as we will discuss in Section [H 

At tree level, the on-shell recursion relations rely on general properties of com- 
plex functions as well as on factorization properties of scattering amplitudes. The 
proof © of the tree-level relations employs a parameter-dependent complex con- 
tinuation u \j,l) n , or "shift", of two of the external massless spinors, j and /, in 
an ?i-point process, 

[j, I) : Xj -> Xj - zXi , Xi^ Xi + zXj . (fO) 

where z is a complex number. The corresponding momenta are then continued 
in the complex plane as well, whereby they remain massless, kj(z) = = kf(z), 
and overall momentum conservation is maintained. 

An on-shell amplitude containing the momenta kj and ki then also becomes 
parameter-dependent, A(z). The physical amplitude is given by ^4(0). When A 
is a tree amplitude or finite one-loop amplitude, A(z) is a rational function of z. 
At tree level, A(z) only has simple poles. These poles arise only from the shifted 
propagators of the amplitude. For example, 

i i 

K r...l.., ^ Kl„ s + z{j-\£r..,\n ' 

if the set of legs {r, . . . , s] contains leg /, which is shifted according to eq. (|10p . 
In the vicinity of the location of the pole z rs , the complex continued amplitude 
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is then schematically given by, 

where h = ±1 labels the helicity of the intermediate state, and the labels L and 
R denote amplitudes with fewer legs, which the propagator eq. (fTTj) connects. 
The number of poles z rs in the complex plane is given by the number of ways the 
set of external legs can be partitioned such that the legs j and I always appear 
on opposite sides of the z-dependent propagator. 
We can now use Cauchy's theorem, 

55 /cT^- 0, (13> 

where the contour C is taken around the circle at infinity, and the integral vanishes 
if the complex continued amplitude A{z) vanishes as z — > oo. Evaluating the 
integral as a sum of residues, we can then solve for the physical amplitude A(0) 
to obtain, 

mo) = - £ b« ^ = = = • ( 14 ) 

poles a r,s h "r -s 

The on-shell amplitudes with fewer legs, Ax and Ar, are evaluated in kinematics 
that have been shifted by eq. (fTUj) with z = z rs , where eq. (fTTj) has a pole, 



K 2 

Zrs = ~{r\£. s ,\n' (15) 

In the following, such shifted, on-shell momenta will be denoted by k(z = z rs ) = 
k. A typical contribution to the sums in eq. (|14p is illustrated in Fig. [TJ 

We have thus succeeded in expressing the n-point amplitude A in terms of 
sums over on-shell, but complex continued, amplitudes with fewer legs, which 
are connected by scalar propagators. These recursion relations can be extended 
to massive QCD and other theories (|46LI47|) . Moreover, for certain helicity con- 
figurations, this recursion relation can be solved explicitly, leading to new all- 
multiplicity expressions for these amplitudes fi8]). 

The basic ingredients to obtain such a recursion relation are complex momenta 
and analysis, which are necessary to make 3-point vertices non-vanishing; factor- 
izability, which is responsible for the simple pole structure; and the vanishing of 
the boundary contribution as z — » oo. At tree level in QCD, one can always find 
complex continuations where this boundary condition vanishes. However, as we 
will see in Section U] below, this is not the case at the one-loop level, and a recur- 
sive approach becomes considerably more complicated. Other theories such as a 
scalar (j) 4 theory have non-vanishing z — > oo behavior already at the tree level, 
which spoils the recursive approach. Studies of the origin of these boundary terms 
and their relation to the Lagrangian can be found in Refs. (|49h . Furthermore, 
additional amplitude structures related to on-shell recursion relations and twistor 
space have been uncovered (f50|) . 
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2.3 Structure of One-Loop Amplitudes 



We now turn to the discussion of gauge-theory one- loop amplitudes, the main 
subject of this review. Prom here on, we will denote one-loop amplitudes explicitly 
with a superscript ^4 1_loo P ) an d tree amplitudes without superscript simply by A. 

Using reduction techniques (|5iy52ll53j) any m-point scalar integral, m > 4, can 
be reduced to scalar integrals with at most four propagators. That is, any one- 
loop n-point amplitude A^ loop can be decomposed into a basis B4 of scalar box, 
triangle, bubble and tadpole integrals, with rational coefficients in four dimen- 
sions. In D dimensions, for example when working in dimensional regularization 
where D = 4— 2e, then the basis Bp is extended to include a scalar pentagon. The 
coefficients of the L>-dimensional basis scalar integrals can be decomposed into 
purely four-dimensional coefficients after expanding in e. Purely rational terms 
are generated when terms higher order in e in the coefficients are multiplied by 
the poles in the integrals, 

4r loop = £ cfif = y: cf= 4 if + n n . (i6) 

Illustrative examples of integrals of the basis B are shown in Fig. [2j 

We will see in the next two sections how to obtain these coefficients and ra- 
tional terms in efficient ways. The scalar integrals contain infrared and ultravi- 
olet divergences that are regulated via dimensional regularization, and depend 
logarithmically or polylogarithmically on momentum-invariants. The integrals 
appearing in eq. (fl~6l) are known and tabulated for example in Refs. (|4"l l54"|) . In 
order to compute one-loop matrix elements the task is therefore reduced to the 
determination of the coefficients, it is not necessary to perform any integrals. 

For amplitudes with only massless particles, the tadpoles vanish. In J\f = 4 
super symmetric theories, as counting of powers of loop momenta in one-loop 
integrals reveals, only box integrals contribute with (D = 4)-dimensional coeffi- 
cients (i.e. free of e terms), and in J\f = 1 supersymmetry bubble, triangle, and 
box integrals contribute, with four-dimensional coefficients @. That is, theories 
with unbroken supersymmetries do not contain purely rational terms that are 
not associated with any of the integrals in the basis. One-loop supersymmetric 
amplitudes can therefore be completely reconstructed from unitarity cuts, as we 
will now discuss. 

3 Extraction of Integral Coefficients via Unitarity 

The problem of computing eq. (|16p has been reduced to determining the coef- 
ficients, Cj, multiplying the known basis integral functions, in the most efficient 
manner possible. The nature of eq. (fTBT) suggests the use of unitarity cuts to 
isolate particular integral coefficients. At the most basic level a unitary cut ef- 
fectively replaces a propagator with an on-shell delta function, i.e. we "cut" the 
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propagator with the replacement, 

1 



p2 _ m 2 _|_ ^ e 



£(+)(p 2 - m 2 ). (17) 



In the original unitarity approach only two propagators were cut but more 
recent developments have highlighted the benefits of applying multiple cuts ([5l 
[12]) . The application of multiple cuts is known as Generalized Unitarity and has 
become the foundation of the most recent developments in the literature. 

Our starting point is the form of eq. (fT6|) decomposed in the basis, we post- 
pone the discussion of the general D-dimensional case using Bry to Section [331 
The purely rational terms are independent of any possible cuts in four dimensions 
and therefore only the remaining "cut-constructible" pieces are accessible via a 
unitarity technique. 

We apply a number of cuts to the expression of the one-loop amplitude and 
match this expression to that of the basis decomposition, eq. (I16p . with the same 
set of cuts applied. This allows us to directly relate the cut expression to the 
basis integral coefficients. This procedure is repeated with as many different sets 
of cuts as is needed to compute all the basis coefficients. Rather than actually 
applying the cuts to the full expression for the one loop amplitude, computed 
for example with Feynman diagrams, we construct the cut expression simply by 
multiplying appropriate on-shell tree amplitudes together. This allows us to take 
advantage of efficient, compact forms of tree amplitudes produced via recursion 
relations, for example those of Section 12.21 

Below, we describe two unitarity approaches for the extraction of integral co- 
efficients. They both utilize knowledge of the integrand, A(l), of the one-loop 
amplitude, A^ loop = J dlA n (l), to derive the basis integral coefficients. The first, 
described in Section 13.11 is based upon the examination of the behavior of the 
loop integrand and the loop momenta in the complex plane. The second, de- 
scribed in Section 13. 2\ known as the Ossola, Papadopoulos and Pittau (OPP) 
method, relies upon computing the coefficients of the general structure of the 
loop integrand itself. 

3.1 Generalized Unitarity in Four Dimensions 

In general we want to isolate as few basis integral coefficients as possible with each 
cut that we consider. It is easy to see that a quadruple cut can be used to isolate, 
on the basis side, a single box coefficient. There are not enough propagators in 
the bubble and triangle integrals to accommodate so many cuts. The set of cuts 
we require to compute all box coefficients corresponds simply to all possible boxes 
that could be present. 

3.1.1 Boxes A very straightforward way to extract a specific box coefficient 
from a quadruple cut expression was proposed by Britto et al. in (|12p . The 
momentum circulating inside a loop without cuts is off-shell and can therefore 



10 Multi-Parton Scattering Amplitudes via On-Shell Methods 

be parametrized in terms of four free components. Applying a cut to one of the 
propagators in the expression of the one-loop amplitude reduces the number of 
free components in this loop momentum by one. Second, third and fourth cuts 
will then reduce the number of free components to zero. The loop momentum 
for the box is then completely frozen by the four delta-function constraints. The 
desired scalar box coefficient can now be read off from the resulting rational 
expression, 



dH d(l) 



J?+ie-KS+(«?) 1 solution, 



(18) 

The coefficient is a product of four tree amplitudes that sit at the four corners 
of the cut box as illustrated in Fig. The momenta flowing into the trees 
satisfy the four cut constraints. In general there are two possible solutions to the 
constraints parameterizing the box. For example with at least one massless leg 
(leg 1) we have (fT0ll55|) . 

(1T|^ 2 ^ 3 ^|1±) 

L± 2<l=F|4t 2 4t 4 |l±> ' liyj 

Inserting each solution into the product of four trees at each corner and then 
summing the two results gives the complete box coefficient, 

rf o = \ E <k. d - = MQMQMQAiiia). (20) 

1 a=± 

Further discussion on the use of both solutions and alternative approaches can 
be found in (56). 



3.1.2 Triangles To compute triangle coefficients, we first apply a triple 
cut to our one-loop amplitude. Unlike the box case, this cut isolates not just a 
single triangle but also any boxes which also contain the same triple cut. Fur- 
thermore we are left with a loop integral with a single free component. Our cut 
expression therefore still depends upon the loop integration. To extract the tri- 
angle coefficient isolated by the triple cut therefore requires two steps. First, we 
must remove the boxes polluting the triple cut expression. Then we must relate 
the triangle integral, which depends upon the remaining free loop-parameter, to 
the scalar triangle basis integral. 

As was first proposed in (|15p . both issues can be solved simultaneously by 
examining the analytic behavior of the triple cut expression in the free parameter 
t of the loop momentum. We choose the following specific parametrization of the 
loop momentum, 

Z" = KY + KY + ^(KYmKY) + ^'-|7 M I^'~>. (21) 
The massless momenta are given by 
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7± = -{K X ■ K 2 ) ± VA, 



det(iQ • A» = (iCi • K 2 f - KlKl 



a = 



S2(gi-7) 
(S1S2-7 2 ) 




(22) 



with 5i = iff. 

By definition, any box terms containing our chosen triple cut will also contain 
an additional propagator, 1/(1 — K4) 2 . Inserting the parametrization for P of 
eq. (|2ip into this additional propagator leads to the development of two poles, 
t±, in the propagator (Z — oc (l/t)(t — t-)(t — t + ). We can take advantage of 
the occurrence of these poles to separate the box terms from the triangle pieces. 

The numerators of these box poles are given by an effective quadruple cut gen- 
erated by the extra pole along with the original triple-cut. As we have seen above, 
the quadruple cut of a box actually corresponds to one of the two contributions in 
the construction of a box coefficient. Therefore each box term can be written as a 
sum of two residue terms, J2i=± di/(xi(t — U)), with di the residue corresponding 
to the box coefficient of the pole U and Xi a constant factor depending upon the 
box in question. 

Analytically, in order to remove these box terms we simply expand the para- 
metrized triple cut integrand expression around t — ► 00. The box terms behave 
as 1/t — > in this limit and thus drop out. Taking a parameter to infinity numer- 
ically is problematic, so instead we use a different approach in the computer code 
BlackHat (10). Considering t as a complex parameter, the box terms appear as 
poles in the complex plane of t, whereas the triangle coefficient is at the origin 
of the t plane. To remove the box terms we simply systematically "clean" the 
complex plane by subtracting all box pole terms from our triple cut expression. 
A final complication to both numerical and analytic approaches is the presence 
of the 1/t factor in the box propagator. To account for this we add back the sum 
of all box terms evaluated at t = , i.e. add (gL — d + )/(xi(t- ~ *+)) f° r each 
box term. An alternative approach to this last step and further discussion on 
the analytic properties of the three- mass triangle can be found in (|15p and ()57j) . 
In addition, the application of on-shell recursion to the computation of certain 
triangle coefficients can be found in (i58j) . 

After the elimination of box terms we are left with a finite power series in t, 



where the Cj are rational coefficients. The upper and lower limits ±n of the sum 
is determined by the theory in question, for example n = 3 for renormalizable 
theories. To relate this sum of terms to the scalar triangle integral we first note 
that, for this particular parametrization, the integrals over any power of t in 
eq. (|23p vanish as proved in (I15p . Then the only remaining term, cq J dt, is in the 
form of a rational coefficient multiplying a scalar triangle integral. In an analytic 
formalism the series expansion around t — > 00 will automatically isolate this term. 




(23) 
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So we can directly relate this sole remaining term to the desired coefficient of the 
basis triangle integral we have isolated with our triple cut. Numerically, our final 
step is to note that, since the power series in t terminates at a finite power, n, 
the full contour integral is equivalent to a discrete Fourier projection (jlOp with 
2n + 1 evaluation points. The triangle coefficient is therefore given by, 

00 = ^ £ C(t e^(^). (24) 

j=-n 

The arbitrary complex number to is the radius of the numerical Fourier projection, 
as illustrated in Fig. [H For technical details we refer to Ref. (|l(jp . 

3.1.3 Bubbles The computation of bubble coefficients proceeds along sim- 
ilar lines as above. A two-particle cut isolates a single bubble coefficient, but will 
also capture triangles and boxes which share the same cut. Again we use the 
differing analytic properties of the terms with additional propagators to separate 
the triangle and box terms from the bubble coefficient. 

As before, we start from a specially chosen parametrization of the loop mo- 
mentum. A two-particle cut leaves two free parameters, y and z. We then choose 
to parametrize the two-particle cut, bubble, loop-momenta as, 

If (y, z) = 1 + (y - I) (i^ - x") + | + ^^(x-^^T) . 

(25) 

Here X: an arbitrary massless momentum, is used to define the massless momen- 
tum = — x^, with its normalization chosen so that K\ ■ x = K 2 /2. 

The triangle and box terms sharing the two-particle cut will contain at least 
one additional propagator 1/(1 — K2) 2 . The parametrization of eq. (|25p will then 
introduce poles in y or z. The residue of the poles in terms of y is given by, 
effectively, a triple cut expression, which is a combination of the pole and the 
original two-particle cut. Contained inside this expression are a single triangle, 
isolated by the effective triple cut, and possibly box terms with the same cut. 

Unfortunately, the more complicated structure of the momentum parametriza- 
tion and the bubble integrand means that a simple extension of the triangle 
procedure is not so straightforward. This is because for the parametrization (|25p 
the integrals over positive powers of y are non-zero, and given by, 

J '*»" = iTTT I (26) 

The integrals over powers of z still vanish. We must therefore alter our approach 
for extracting the bubble coefficient. As before, we discard any terms in the ex- 
pansion around z — > 00 that depend upon z, but we must retain the coefficients of 
all powers of y (which are all guaranteed to be positive due to our parametrization 
choice). For a renormalizable theory the maximum power is y 2 . We can relate 
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these terms to the scalar integral using eq. (|26p . which is independent of y and z, 
J dydz. The sum of the resulting terms forms only part of the bubble coefficient. 

The source of the remaining contribution is the triangle expression from the 
residues of the poles in z. This effective triple cut does not vanish in the double 
expansion in y and z. Therefore, we obtain an additional contribution which 
needs to be subtracted from the part of the bubble coefficient computed in the 
previous paragraph. Further details on this subtraction term can be found in 
Ref. (fT5|) . 

To numerically extract the bubble coefficient we must first "clean" the complex 
plane of all pole terms. Computing the residues of each pole involves simply 
computing the triple cut at the location of each pole. Once all pole terms have 
been removed we are free to extract the bubble coefficient from the remaining non- 
pole terms using a double discrete Fourier projection, in both y and z. Instead 
of naively evaluating at as many points as there are coefficients in y and z we use 
the nature of our parametrization to reduce the number of points at which we 
need to evaluate z by one. The coefficient in a renormalizable theory is given by 



-t 

20 ^ 

3=0 



B (y = 0, z = t e 2nij/5 ) + 3B (y = 2/3, z = t e 2nij/5 ) 



(27) 



B(y,z) denotes the two-particle cut from which triangles and boxes have been 
subtracted that share this cut. The more general case is given in ()10p . Alternative 
approaches to the computation of the bubble coefficients have been proposed in 
the literature, such as Mastrolia's use of Stokes Theorem (|20p. 

3.1.4 Massive particles The expressions we have given above are for 
amplitudes with purely massless particles. The addition of massive particles 
which do not circulate in the loop is straightforwardly accommodated within the 
above, with no changes. Including massive particles inside the circulating loop 
requires further exposition. Two changes are required, firstly the loop momentum 
parametrization needs to be extended to include massive particles. Secondly we 
must also compute the coefficients of tadpole integrals in addition to the bubbles, 
triangles and boxes. A detailed discussion of the extension of the procedure to the 
computation of massive particles as well as the computation of the tadpole terms 
themselves is given in (59). In addition there is the problem of the wave-function 
renormalization with massive particles, which has been addressed in Ref. (poli - 



sh Extraction of Integral Coefficients at the Integrand Level 

So far we have described methods which used analytic limits or a combination of 
the subtraction of poles on the complex plane and discrete Fourier projections. 
An alternative approach developed by Ossola, Papadopoulos and Pittau (OPP) 
uses the knowledge of the general structure of the integrand A(l) instead. The 
integrand is built up from a standard set of terms. These terms either vanish 
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upon integration or correspond to one of the scalar integral basis functions of 
eq. (fl~6|) . Computing the scalar basis integral coefficients then reduces to the 
problem of finding the coefficients of the OPP decomposition of the integrand. 
In four dimensions, the integrand can be written as, 

[) ~ ^ D h D i2 D Vi D iA + A1A2A3 

l<n<i 2 <n n i2 l<ii<n 11 

with the propagator Di = (I — K^) 2 — mf, where the mass of the cut propagator 
with momentum /j has mass rrij. From now on we include massive particles in 
the discussion. The form of the numerators in eq. (|28p depend upon the basis we 
choose for the loop momenta. We wish to choose this momentum basis so that 
each scalar integral basis coefficient of B4 corresponds to a single term in the 
integrand decomposition and the integrals over the remaining structures vanish. 
In the form originally presented by OPP (|13p . a momentum parametrization for 
the box, triangle and bubble very similar to eqs. (fl~9j). (f2T|) and (f25j) was used. 

An alternative momentum parametrization, presented by Ellis et al. (|6ip . is 
related to the van Neerven-Vermaseren basis (|52p . The generic form of a momen- 
tum in this basis is 

Dp D T 

^ = E< + E«^- ( 29 ) 

3=1 i=i 

This is a decomposition into to two sets of basis vectors. The vectors span the 
physical space defined by the external legs Ki and the vectors span the space 
transverse to this physical space. For a box Dj- = 1 and Dp = 3, for a triangle 
At = 2 and Dp = 2 and for a bubble Dt = 3 and Dp = 1. The basis vectors 
are chosen such that i%i ■ nj = 5ij, nj • K j = and Ui-Vj = 0. The Vj's are chosen 
such that any cut legs are on-shell. 

The numerators of the propagator terms are then arranged in the following 
way (f6T|) . 

^1112*3*4(0 = C *l*2*3*4 + ^1*213*4*1' (^0) 
^11213(0 = C *l*2*3 + C *l*2*3*l + C *l*2*3* 2 + C *l*2*3(*l ~~ ^2) (31) 
+^1^2(Cj 1 j 2 j 3 + Ci^i^l + c i 1 i 2 j3^2) 5 

WO = 4*2 +4*2*1 +4*2*2 + 4*2*3 + 4*2 (*?-*!) (32) 

+4*2 (*2 - *l) + 4*2*1*2 + 4*2*1*3 + 4*2*2*3- 

MO = 4 +4*1 + 4*2 + 4*3 + 4*4 > (33) 



for the tadpole, bubble, triangle and box, i = 1,2,3,4 respectively. Here the 
tj = (nj ■ I) depend on the specific parametrization of the loop momentum, which 
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is chosen differently for boxes, triangles, bubbles and tadpoles. The coefficients we 
wish to compute are the c° terms which correspond to the scalar basis integral 
coefficients. We will not discuss the computation of the tadpole terms here, 
methods to compute these can be found in p3l[T9ll59ll60l[6T]) . The remaining 
coefficients multiply tj terms (or combinations of such terms), which vanish upon 
integration over I. As in section \3A\ the choice of the representation for the loop 
momentum is crucial in this approach. 

Trying to solve for the entire set of coefficients at once is clearly not the optimal 
approach. Using unitarity cuts we can isolate individual terms of the integrand. 
The application of all possible cuts allows us to systematically solve for all of the 
coefficients sequentially. Isolating a single box term with a quadruple cut 

works in the same way as in section 13.11 Again we have two solutions for the 
completely frozen box loop momentum. For the momentum basis corresponding 
to the integrand decomposition eq. (i30j) these are 

1* = v ? + v% + vg ± ^Jm 2 - ( Vl +v 2 + v 3 ) 2 <. (34) 

As above, the box coefficients in the OPP basis are given by products of four tree 
amplitudes, 

c ° , 2 i 3 i4 = \ E d ^ d * = Mla)A 2 (l a )A 3 (l a )A 4 (l a ) , (35) 

Z a=± 

c hi2hu = 2 ( d+ ■ 

Here, c\ li2i ^ iA is the desired scalar box integral coefficient, do (cf. eq. (f20j) ). The 
other coefficient is needed for the computation of the triangle coefficients. 

For the computation of the triangle coefficients we apply a triple cut to isolate 
the particular triangle we are interested in. Much like before this means that we 
also have box terms polluting the result. To solve this problem we subtract the 
complete box contribution from the integrand. So in order to find the triangle 
coefficients cf lia j 3 we do not evaluate A\(l)A2(l)A 3 (l), but instead, 

A 1 (l)A 2 (l)A 3 (l)-Y / dlll t 3uil \ (36) 

at different choices for the loop momentum of the triangle. This is in contrast to 
the procedure of the previous section. There we subtracted the residue of the box 
pole from the triangle rather than the entire contribution of the box (]20p from 
the integrand as above. The triangle momentum parametrization corresponding 
to the integrand structure eq. (|3ip contains one free parameter which we label a, 

P = vl + v% + y/ml -a 2 - (ui + v 2 ) 2 < + a n£. (37) 

There are seven unknown coefficients in eq. (|31[) . so we need to evaluate this 
at seven, in principle arbitrary, values of a to determine all coefficients. For 
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increased numerical stability, however, they can be chosen to lie on a circle as in 
the previous section. The resulting set of linear equations can then be solved to 
find the full set of coefficients, which is needed in the computation of the bubble 
coefficient. 

Finally, in order to compute the bubbles we apply a double cut, which again 
isolates a single bubble but also any triangle and box coefficients that share the 
same cut. Again we simply subtract the unwanted terms from the integrand to 
remove them. To find the coefficients ,■ we evaluate 



at different values of the loop momentum of the bubble. Here the bubble loop 
momentum corresponding to the parametrization eq. (|32p is, 



with two free parameters, ot\ and a.2- Solving for the nine coefficients in eq. (132j) 
requires nine linear equations. These can be generated by choosing nine different 
values of I via the choice of the two free parameters. The full set of nine coefficients 
is only required if we wish to compute the tadpole coefficients, which only appear 
if massive particles are present in the loop. Further details can be found in 
Refs. (JSaEU). 

3.3 Rational Terms from D-Dimensional Generalized Unitarity 

We have so far restricted ourselves to working with cut legs in four dimensions, 
keeping the rational terms out of reach. Considering cuts in D dimensions allows 
us to use the Bp integral basis of eq. (|16p. bringing the computation of all terms 
within our grasp. This is related to van Neerven's important observation that 
dispersion relations for Feynman integrals converge in dimensional regulariza- 
tion (|62p . There are multiple different approaches to go beyond four-dimensional 
cuts as already proposed in Refs. (j22"j) . The additional D — 4 components can 
be related to massive terms in four dimensions ()22j) . In a related approach the 
D — 4 components can be converted to an additional integral which can be used 
to compute the full amplitude (I63h . 

Alternatively, we can directly extend the approaches of either Section 13.11 or 
Section 13.21 as we will now describe, starting with the generalization of Sec- 
tion E21 

3.3.1 D-DlMENSIONAL UNITARITY AT THE INTEGRAND LEVEL OPP have 
suggested a two-step computational procedure. Here, the D — 4 terms in the nu- 
merator of the integrand are computed using a separate set of Feynman diagrams. 
The corresponding Feynman rules have been derived for QCD and electroweak 




(38) 




(39) 
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theories in a series of papers (|24^I64|) . The D — 4 contributions from the denomi- 
nator of the integrand are found by computing coefficients of an extended OPP 
basis structure for the integrand (|13 y24p . 

A direct extension of the OPP approach to D dimensions was derived by Giele 
et al. (21). This approach combines trees in higher dimensions with an extended 
integrand basis for the one-loop integrand, taking into account the additional 
structure in higher dimensions. Again as in Section 13.21 the choice of the loop 
momentum determines the form of the integrand structures. The form of the 
higher dimensional loop momenta can be decomposed into a four-dimensional 
part I and an orthogonal (D — 4)-dimensional part I, P = 1^ + 1^. The on-shell 
constraint then means that I 2 = I 2 + I 2 and so the four-dimensional component 
effectively becomes massive with mass —I 2 = fi 2 , where u 2 is the "scale" of the 
higher-dimensional subspace. 

Since the external momenta remain in four dimensions, we find that there are a 
limited number of additional integrand structures that can be present because the 
external momenta are orthogonal to the D — 4 additional transverse dimensions. 
The numerator structures can only depend upon the higher-dimensional terms 
through even powers of fi. Similarly, there can be at most one additional cut leg, 
since fi 2 is fixed by the fifth constraint. Therefore we must include a pentagon 
in the integrand basis and so for a renormalizable theory the extended integrand 
structures are given by, 

*n{i) - 2^ TTTTTr~n~n~ + 2^ 



*ll«2»3«4 

I<ii<i2<«3<«4<i5<n ~*l~*2~»3~*4~»5 l<i t <i 2 <i 3 <i 4 <n 

_l_ ^nnhW _|_ bj r i 2 {l) ^ ^ 5^ (I) 

l<h<i 2 <ia<n D h D i2 D ki l<i 1<i2 < n D h D i2 l< ix < n D h 

Here the numerator coefficients are given by, 



-ll«2»3*4«5 
d, 



H«2«3M V 
^«l«2i3 ( 
^ill2 ( 
5*1 ( 



c° • • ■ ■ (41) 

°Jl22«3*425 ' \ J 

^1121314(0 + A* ( c hi2i3i4 h c iii2hi<t) M c iii2i3i4' 
c il 1213(0 + A* (^1 c u«2«3 ^ 2 c ixi 2 i 3 + c hi 2 

6i lia (l) + ii 2 c}% 2 , 
Oil (0 • 



where again tj = rij ■ I. 

The attentive reader may be worried about how this basis can be used if we 
have a fractional number of additional dimensions, for example when D = 4 — 2e. 
This problem can be sidestepped by noting that the dependence of the amplitudes 
upon the {D — 4)-dimensional subspace is linear. So the fractional dimensional 
structure of the amplitude can be reconstructed by evaluating the amplitude at 
two different integer dimensions. Using A^~ loop = Aq~ 1oop + (D — 4)^4}~ loop the 
full D-dimensional amplitude can then be reconstructed (|2ip . This allows us to 
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freely choose the exact form of the regularization scheme. The integer dimension 
chosen when evaluating the trees that enter the calculations is arbitrary up to the 
constraint that the dimension has to be an even number if fermions are present 
in the loop. 

The final step is to relate the new integrand structures of eq. (I4ip , proportional 
to powers of fi, to their contribution in eq. (|16p . As before, terms proportional to 
tj terms will vanish and so we are left with terms proportional only to fi. Scalar 
integrals multiplied by powers of fi can be related to scalar integrals in higher 
dimensions. From eq. (|41|) we have the following different integrals, 

dPi ^ = if+L.^o 



-Di 1 Di 2 Di ll Di4 2 



H1213M 



,p, ^ (£>-4)(£>-2) D+4 D ^_l 

4 6 

d Dj M 2 ( D -^) rD+2 D^i_l 

D h D i2 D i3 2 ^ 2 

jDj V? ( D ~ 4 ) tD+2 K + m l 1 q Uo , 

d D h D i2 - 2 '™ 2 6 bl1 ' 12 ' {4Z) 

where is the mass of the bubble labelled by i\ and %2- Examining the limits 
of these integrals as we return to four dimensions, we see that the integrals are 
finite and correspond to purely rational numbers. Each new integrand factor 
therefore contributes to the rational term. 

Computing the one-loop amplitude is now very similar to the cut-computation 
procedure of Section I3.2L Starting with the pentagon we find that the loop 
momenta will be completely frozen by the four cuts on the four-dimensional 
component of l§ and the constraint on /j? = —I 2 from the fifth propagator (7 + 
I — -K5) 2 = m|. The resulting expression is then inserted into the penta-cut so 
that we find the pentagon coefficient using, 

c° ili2i3Ui5 = A!{l)A 2 {l)Ml)Ml)Ml)- (43) 

The trees must be evaluated, as described above, at two different integer di- 
mension choices for the internal loop momentum. The pentagon coefficient is 
proportional to D — 4 and so will vanish in the four-dimensional limit. We need 
it only for subtraction when computing the box terms as we will now explain. 
The box is computed as before but with the pentagon subtracted, 

^(0^2(0^3(0^4(0 - E etlt T ts(0 - (44) 

With three additional coefficients to be determined, compared to the four-dimen- 
sional case, we need three additional evaluation points in order to generate enough 
equations to solve for all coefficients. All coefficients are required for the com- 
putation of the triangle terms. Computing the triangle and bubble contributions 
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to the rational terms follows in a similar vein as for the boxes and we direct the 
reader to Ref. (|2ip for further details on the computation of all such contributions. 

3.3.2 Four-dimensional Generalized Unitarity with a D — 4-Dimen- 
SIONAL Mass The differing dependence of the rational terms of eq. (|41|) on 
u 2 suggests that we could use the analytic structure of the integrand to compute 
the rational terms at the integral level. This approach of Badger (23) is similar 
in spirit to the extraction of the cut-coefficients described in Section 13. 11 Here 
the loop momentum is massive and four-dimensional rather than massless and 
.D-dimensional . 

To compute the rational contribution of the box, we start from the quadruple 
cut expression of a box. Now each cut leg has an additional mass n 2 . This cut 
isolates a single box coefficient and also any pentagon terms which share the same 
cut. Schematically we have, 

^ 1 (Z(M 2 ))A 2 (/(^))A 3 (/(^))^ 4 (/(/i 2 )) = n + r 2 n. 2 + + 7 f_ 2 V (45) 

j Xi fa ) 

where /i 2 is the pole in u 2 for the pentagon i, \i is a constant factor depending 
on the pentagon in question. The additional propagator of the pentagon terms 
shows up as a 1 /(/j, 2 — /i 2 ) factor. We wish to separate these terms as well as 
the other terms in the powers series in u 2 from each other. Here we only want 
the coefficient r% of /x 4 , which is the term that contributes to the purely rational 
part, cf. eq. (g2]). So simply expanding A 1 {l{n 2 ))A 2 {l{^ 2 ))A z (l(ii 2 ))A A {l{^ 2 )) / ^ 
around the limit n 2 — > oo will give us the coefficient directly. 

Moving on, performing a triple cut isolates a single triangle as well as boxes 
and pentagons. These objects contain poles in both n 2 and the unconstrained 
component of the loop momenta, t. We are only interested in the coefficient of 
/j, 2 and so series expanding around both t — > oo and u 2 — > oo isolates this single 
term. For the bubbles a similar procedure applies now in three parameters, a 
complete description is given in (|23|) . 

As in Section [3"7Ll taking an infinite limit numerically is difficult to do. Similar to 
the numerical approach to the computation of the cut-containing terms described 
in Sect. 13.11 we can adapt Badger's method. This is the approach implemented 
within BlackHat (|10p . Starting from eq. (|45p we "clean" the complex [i plane by 
subtracting all pentagons at their poles. We then perform a discrete Fourier pro- 
jection in |U to extract the rational contribution r$. For the rational contribution 
of a triangle we consider the triple cut expression and then clean the complex /i 
plane by subtracting all terms with an additional propagator. Each such term 
has a pole in t and its numerator is a quadruple cut. Here we are subtracting the 
quadruple cut residue of the additional propagator and not the integrand box and 
pentagon terms, as we would in the OPP approaches. We sidestep therefore any 
loss of numerical stability that would arise from cancellations between these box 
and pentagon terms since these pieces are never separated out in the quadruple 
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cut. The computation of the rational contribution from the two-particle cut of a 
bubble follows along a similar line. 

4 On-Shell Recursion at One Loop 

As above, an alternative approach to obtain the full amplitude is to use four- 
dimensional unitarity methods for the cut-constructible terms and to use recur- 
sion relations for the purely rational remainder. The recursive approach devel- 
oped in Refs. (J5JED has already been reviewed extensively in ([3]). While this 
approach is very useful for analytical calculations and even allows to obtain 
closed- form all- multiplicity results (|65p . it is not straightforward to cast into 
a form suitable for numerical implementation. This is due to the removal of 
spurious singularities via the introduction of overlap terms, which is difficult to 
perform in an automated fashion. The approach has subsequently been modi- 
fied to allow efficient numerical implementation into the automated C+- 1- library 
BlackHat (fT0|) . We will discuss here this modified approach, for amplitudes with 
massless particles in the loop. 

We begin by dividing the amplitude into cut-constructible and rational terms, 
as in eq. f)16[) . The rational terms are defined by setting all scalar integrals to 
zero, 

4T 1o ° P = E cf= 4 l? + K n , (46) 

C n = E C f =il F > = ^- l0 ° P TD n • (47) 

jeB 4 



In the following we assume that the cut-containing terms have been computed 
via the generalized unitarity methods described in the previous section. 

The basic idea is to complex continue 1Z and use Cauchy's theorem to recon- 
struct the rational term from its poles in the complex plane, similarly to the 
tree-level approach introduced in Section I2.21 We add this rational term to the 
previously-computed cut terms C to obtain the full physical amplitude at z = 0. 
However, two of the basic premises of the derivation in Section I2.2I do not hold 
in general at one loop: For one, we cannot always find shifts (|10p such that the 
amplitude vanishes as z — > oo. And for another, the division of the amplitude 
into cut and rational parts introduces the presence of spurious, unphysical poles 
in the complex plane when considering the rational part separately. That is, the 
rational part has the following form upon complex continuation, 

n n ( z ) = n p n (z) + n s n {z) + n^ z , (48) 



where 7£„ denotes the contribution from physical poles, which are as at the tree 
level simple poles, whereas the spurious poles can be simple or double poles, and 
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the contributions TZ^' S ° Z display polynomial growth in z: 

= E^V> (49) 

= E fr^S + . 

y^(z-^) 2 2-2/3; 

Cmax 

cr=0 

Here, ^4, -B, C, and D are functions of the external momenta. We postpone the 
case where TZ^(z) has a more complicated structure to Section [4,21 

The physical contribution is recursively computed as in the tree level case, 

E = EE(^ = ^)^-^ = ^ (50) 

poles a " r,s h I r - s 

+ = Z rs )—^A^ h ( Z = Z rs ) 

"-T...S 

+ A h L {z = Z rs ) l ^.- s) A~ R h {z = Zrs)} ■ 

Here, we have expressed the rational term as a sum of products of rational terms 
from lower-point amplitudes (defined according to eq. (|47|) ) with lower-point tree 
amplitudes. The last term with T corresponds to a one-loop correction to the 
propagator. Eq. ([50]) is illustrated in Fig. [SJ 

The full amplitude is found by combining all rational and cut contributions, 

Ai- lo °P(0) = C n (0) - y Res ^^l + n^(0)+n^ ez . (51) 

poles a 

We will now show how to obtain the last two contributions, starting with the 
spurious pole contribution. 

4.1 Spurious Poles 

The division into cut and rational parts introduces spurious singularities in each 
of these terms which however cancel in the full amplitude. These spurious sin- 
gularities are already present in real kinematics. Cauchy's theorem requires us 
to sum over all poles, whether physical or unphysical. In Ref. ([9]) this problem 
was remedied by adding additional rational terms to the cut part. These rational 
terms are constructed such that the cut and the rational parts individually do not 
contain spurious singularities upon continuation into the complex plane. How- 
ever, these extra rational terms then contribute additional terms to the residues 
at the physical poles in eq. (]50p which have to be subtracted from the rational 
part 1Z in the recursive construction by so-called overlap terms in order to avoid 
double counting ([9]). This approach leads to compact expressions in analytical 
calculations but is not particularly amenable to numerical implementation. 
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An alternative way of dealing with the spurious singularities is to make use of 
the fact that we know that they cancel in the full amplitude. In other words, we 
can extract the spurious residues from the known cut parts, 

k» s (°) = - E n^ = + £ 25 nr- (52) 

spur poles f3 spur poles /3 

where C n (z) is the complex continued cut part. The spurious poles in C n (z) come 
from the vanishing of complex continued Gram determinants, A(z), associated 
with bubble, triangle, and box integrals. Since the spurious poles cancel between 
the rational and the cut parts, the spurious contribution to the residues from the 
cut part can only be rational. To compute the residue we therefore series expand 
the logarithms and polylogarithms around the location of the vanishing Gram 
determinants and obtain a series of rational functions. The spurious contribution 
is thus given by evaluating 



7Z*(0) = E J* 

A m ( 2:/3 )=0 



c? =i < 



E — > ( 53 ) 



A m ( Z(3 )=0 



where the subscript "rat" indicates that we take the rational part of the series 
expansion of the integrals around the spurious poles. We have introduced the 
abbreviation E@ for these rational terms. The spurious poles zp are located where 
the shifted Gram determinants vanish, A m (zp) = 0, with m = 2, 3,4 for bubble, 
triangle, and box integrals, respectively. Note that poles in Gram determinants 
of box integrals will in general also appear in the daughter triangle and bubble 
integral^. Singularities in triangle integrals will feed down into the daughter 
bubble coefficients, but not affect the parent box integrals. 

The expansion of an integral around the location where its Gram determinant 
vanishes can be obtained (|66p by using the "dimension-shifting relation" (j53j) 
iteratively, 

X F = \ ( £ c * J f-i W + <j-l-D)co If + 2 J , (54) 

where denotes the lower-point integral obtained from Xf by removing the 

ith propagator. The coefficients Cj and cq are given by 

3 3 

Ci = E(^ _1 ). fe , c = J2c i ~A j , (55) 

fe=l k=l 

where Y -1 denotes the inverse of the modified Cayley matrix. These modified 
Cayley matrices are listed for example in (1541) . As indicated, the coefficient cq is 
proportional to the Gram determinant. We obtain the series expansion of if in 



3 By daughter integrals we mean integrals that can be found from the parent integral by 
collapsing one or more of the loop propagators. 
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terms of its Gram determinant, 



1 D 



rat 



r*A*, (56) 

A;=0 

with rational coefficients that are found by using eq. (|54p iteratively. The 
explicit expression for the rational expansion of the three-mass triangle is listed 
in Ref. ()10j) . Similar expansions for the remaining integrals can be obtained as 
described above and will be listed in a forthcoming publication (C.F. Berger et 
al., in preparation). 

Numerically, we can evaluate eq. (|53p by using a discrete Fourier sum, which 
here only approximates the contour integral. We evaluate the quantity given in 
the square bracket of ([53]) . E@(z)/z, at m points equally spaced around a circle 
of radius 5p in the z plane, centered on the pole location zp, 

*n(0)~-EI> Zp + Sp e^ ■ (57) 

The sum over (3 runs over all locations of spurious poles where Gram determinants 
vanish. For technical details of the numerical implementation we refer the reader 
to Ref. (HDD. 



4.2 Contribution from Infinity 

The remaining rational contribution TZ^ rs ° z is the boundary contribution in the 
contour integral as z — > oo. Although it is generically possible to find com- 
plex continuations that have vanishing boundary contributions, these shifts have 
in general additional contributions that cannot be constructed recursively as in 
eq. ([50]) . that is, we have instead, 

7^(0) = ^recursive + ^P nonstop ^Procursive _ _ J- ^ KjA . ( 58 ) 

l a Z ~ Z P % 
poles p 

These 'non-standard' contributions (labeled by the superscript "nonstd") arise 
in configurations where two external momenta with the same helicity are on 
one side of the partition and all other legs are on the other side. The complex 
factorization properties of these configurations are not yet fully understood. The 
sum over poles (3 in eq. (J58]) is only over the channels that factorize, i. e. that do 
not display the aforementioned problematic behavior. Conversely, it is possible 
to find shifts that avoid these 'non-standard channels', however, at the price of 
reintroducing a boundary contribution. 

The solution to this problem, developed in Ref. ([9|), is to use two independent 
complex continuations to determine the boundary contribution. Let us denote the 
primary shift by [j, I) and the auxiliary shift by [a, b) in the notation of eq. (|10p . 
We then have two relations for the same amplitude, analogous to eq. (|5T|) . 

^- loop (0) = Inf ^ loop + C n (0) - Inf C n + ^ ^cursive, [7,0( ) j ( 59 ) 
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^l- lo °P(0) = C„(0) - Inf C n + ^recursive, [a,b) + nonstd,[a,fe> + ^5 [a,6) ( Q ) _ (g ) 

[a,b) 

We have indicated with additional superscripts which shift has been employed. 
11% is evaluated according to eq. ([53]) . We have also used that, 



^large z = j nf ^1-loop _ j nf (j /g^ 

W) W) 

where Inf A is the unknown large- z behavior of the full amplitude found from a 

Laurent expansion of A l - lo °v(z) around z = oo with the shift [j, I), and similarly 
for C n (z). Eqs. (|59[) and (|60p thus both contain unknown terms. If we now apply 
the primary shift [j,l] to the auxiliary recursion (|60p . and take the limit 2 — > 00, 
we can extract the large- z behavior of the primary shift, 

Inf A3,- 1oo p = InfC n -Inf ( Inf C n ]+In£ ( Vl rccursivc [a ' b) )+lni ( K^ a ' b Uo la b) 

(62) 

where now all terms on the right-hand side are either known or recursively con- 
structible, if 

Inf ^nonstd,[a,^ =Q ^ 

Putting everything together, we find the full amplitude from eq. (|59p . with 
([50]) . ([53]) . and ([62]) . and the cut terms constructed as described in the previous 
section. 



5 Conclusions and Outlook 

In this review we have presented an overview of recent developments in the calcu- 
lation of multi-parton scattering amplitudes at the one-loop level. These devel- 
opments are based on on-shell techniques that make efficient use of the physical 
properties of the hard scattering, such as unitarity and factorization. The basic 
ingredients in these new approaches are on-shell tree-level or lower-point one- 
loop amplitudes instead of Feynman diagrams, thus sidestepping many of the 
complications associated with the use of Feynman diagrams. 

Furthermore, these new techniques allow for efficient algorithmic implemen- 
tation and hence the construction of efficient, numerically stable, and fast com- 
puter codes, such as BlackHat (QUI), CutTools/OneLOop Rocket ([27]) . 
and others (|28p . With these new techniques and computer tools a flurry of re- 
sults relevant for the LHC has recently been computed (|3i y 32 p 33,34,35,36). and 
we expect further rapid progress in the near future. 

Nevertheless, much work remains to be done to bring one-loop calculations 
to the same level of automatization as tree- level computations, ideally starting 
from a Lagrangian and producing complete events including parton shower and 
hadronization corrections. Further open issues include, for example, a better 
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understanding of the complex factorization properties of one-loop amplitudes and 
the generalization of the new techniques to higher loops in nonsupersymmetric 
theories. 

In summary, the last few years have seen an unprecedented progress in the 
development of techniques for the computation of multi-parton one-loop scatter- 
ing amplitudes which are an essential ingredient in precision calculations for the 
LHC. These new methods have also been used to study the higher- loop structure 
of J\f = 4 supersymmetric Yang-Mills theory and N = 8 supergravity ()67|) . Their 
basic ingredients are unitarity, factorization and complex analysis, properties that 
are quite generic. It is thus not inconceivable that these new techniques will find 
further application beyond those presented or referenced in this review. 
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Figure 1: Schematic representation of a typical tree recursive contribution to 
eq. (|14|) . The labels £ T" refer to on-shell tree amplitudes. The momenta j and I 
are complex continued, on-shell momenta. 




(a) (b) (c) (d) 

Figure 2: Representative examples of integrals appearing in eq. (|16[) : (a) a box 
integral, that is, a 4-point integral, (b) a triangle (3-point) integral, (c) a bub- 
ble (2-point) integral, and (d) a tadpole integral. Each corner can have one or 
more external momenta attached to it. The tadpole vanishes when all internal 
propagators are massless. 




Figure 3: A quadruple cut one-loop integral isolating the single box term 
MKlKlKlKl). 
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Figure 4: The points on the circle used by the discrete Fourier projection, cf. 
eq. dMD- 




The labels 'T' and 'L' refer to tree and the rational part of loop vertices, respec- 
tively. 



